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The calculation of QCD corrections to the width difference AF in the _Bs -meson system is pre- 
CL|| sented. The next-to-leading order corrections reduce the dependence on the renormalization scale 

p • significantly and allow for a meaningful use of hadronic matrix elements from lattice gauge theory. 

^ ■ At present the uncertainty of the lattice calculations limits the prediction of AF. The presented work 

has been performed in collaboration with Martin Beneke, Gerhard Buchalla, Christoph Greub and 

Alexander Lenz. 



1 Introduction 

Precision analyses of flavor-changing transitions are of experimental top priority in the forthcoming 
years. Decays of B mesons provide an especially fertile testing ground for various reasons: they 
allow for a high precision determination of three of the four parameters characterizing the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix [[l|], including the CP- violating phase 7. Since flavor-changing 
transitions of i? mesons are always suppressed by small CKM elements and heavy electroweak gauge 
boson masses, it is well possible that B physics experiments will reveal new physics. The large mass 
rrib of the 6-quark further allows us to control hadronic uncertainties. Fermilab's CDF, DO [0] and 
the planned BTeV [Q experiment prepare a dedicated B physics program. Other studies are in 
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Figure 1: Bg — Bg mixing in the Standard Model. The zigzag lines represent Ty -bosons or charged 
pseudo-goldstone bosons. 



progress at CLEO, LEP and at HERA-B [^ or planned for the future LHCb [Q] experiment. While 
5-factories [Q] only produce Bd^B^ and B^ mesons, LEP and the hadron colliders also provide Bg 
mesons. Like their K, D and B^ counterparts Bg mesons mix with their antiparticles. Therefore 
the two mass eigenstates Bh and B^ (for "heavy" and "light") are linear combinations of Bg and 
Bg and differ in their mass and width. In the Standard Model Bg — Bg mixing is described in the 
lowest order by the box diagram depicted in Fig. [I]. The dispersive part of the Bg — Bg mixing 
amplitude is called M12. In the Standard Model it is dominated by box diagrams with internal top 
quarks. The absorptive part is denoted by ri2 and mainly stems from box diagrams with light charm 
quarks. ri2 is generated by decays into final states which are common to Bg and Bg. While M12 can 
receive sizable corrections from new physics, ri2 is induced by the CKM-favored tree-level decay 
h -^ CCS and insensitive to new physics. Experimentally Bg — Bg mixing manifests itself in damped 
oscillations between the Bg and Bg states which are governed by M12 — iri2/2. We denote the mass 
and width differences between Bh and Bl by 

Am = Mh - Ml , AF = T^ - Th- 

By solving the eigenvalue problem of M12 — iTi2/2 one can relate Am and AF to M12 and F12: 
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Am equals the Bg—Bg oscillation frequency and has not been measured yet, but we know the lower 
bound Am > 14.9 ps^^ from LEP data [^. It can be shown that this bound implies IF12I/IA/12I ^ 



2 QCD effects 



0.01. In deriving (jl]) terms of order |ri2/Mi2p have been neglected. in (|^) is a CP-violating 
phase, which is tiny in the Standard Model, so that ATsm = 2|ri2|. In the presence of new physics 
argMi2 and thereby (p can assume any value, (p can be measured from CP-asymmetries, which 
requires the resolution of the rapid Bg—Bg oscillations and tagging, i.e. the discrimination between 
Bs and Bs mesons at the time t = of their production. From @) one verifies that a non-vanishing 
also affects AF, which can be measured from untagged data samples and therefore involves better 
efficiencies than tagged studies. Unlike in the case of Bd mesons, the Standard Model predicts a 
sizable width difference AF in the Bg system, roughly between 5 and 30% of the average total width 
F = (F^, + F//)/2. Now the decay of an untagged Bg meson into the final state / is in general 
governed by two exponentials: 

F[/,t] oc e-^-'\{f\BL)\' + e'-^"'\{f\Bn)\\ (3) 

If / is a flavor-specific final state like D~7r~^ or Xi'^u, the coefficients of the two exponentials in (|]) 
are equal. A fit of the corresponding decay distribution to a single exponential then determines the 
average width F up to corrections of order (AF)^/F. In the Standard Model CP violation in B^—Bg 
mixing is negligible, so that we can simultaneously choose Bl and Bh to be CP eigenstates and 
the b — > CCS decay to conserve CP. Then B^ is CP-odd and cannot decay into a CP-even double- 
charm final state fcp+ like {J / 1^4)) l=o,2, where L denotes the quantum number of the orbital angular 
momentum. Thus a measurement of the Bg width in Bg -^ fcp+ determines F^. By comparing 
the two measurements one finds AF/2. In the presence of a non-zero CP-violating phase 4> this 
procedure measures [^ 

AFcos0 = AF5mcos^0- (4) 

The extra factor of cos (p stems from the fact that in the presence of CP violation both Bi and Bh 
can decay into fcp+- CDF will perform this measurement with Bg — > Djn^ and Bg — > J/V^^ in 
Run-II of the Tevatron [|^]. 

2 QCD effects 

The Bg—Bg mixing amplitude of Fig. |l]and the Bg decay amplitude are affected by strong interaction 
effects. AF^M = 2|Fi2| involves various different mass scales and the the QCD corrections associ- 
ated with these scales require different treatments. In the first step an operator product expansion at 
the scale Mw is performed to integrate out the W^-boson. The Standard Model b -^ ccs amplitude is 
matched to its counterpart in an effective field theory in which AB = 1 transitions (B is the bottom 
number) are described by four-quark operators. The corresponding effective hamiltonian reads 

K// = ^K;^c. ( E CrQr + CsQs) , (5) 
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with the operators 

Ql = {biCj)v-A{CjSi)v-A Q2 = {biCi)v-A{CjSj)v-A, (6) 

Qs = {piSi)v-A{qjqj)v-A Qi = (PiSj)v-A{qjqi)v-A, (V) 

Qb = {biSi)v-A{qjqj)v+A Qe = {biSj)v~A{qjqi)v+A, (8) 

Q, = £^m,ka^'^{l~j,)^s,G;^. (9) 

Here the i,j are colour indices and a summation over q = u, d, s, c, b is implied. V ± A refers 
to 7'^(1 ± 75) and S — P (which we need below) to (1 — 75). The current-current operators Qi 
and Q2 stem from M/-boson exchange between the be and cs lines. Qs-e are four-quark penguin 
operators and Qg is the chromomagnetic penguin operator. The Wilson coefficients Ci contain the 
short-distance physics and are functions of the heavy W and top quark masses. Since they do not 
depend on long-distance QCD effect, they can be calculated in perturbation theory. C^-e are very 
small. The matching calculation determines the C,;'s at a high renormalization scale /i = 0{Mw)- 
The renormalization group (RG) evolution of the coefficients down to fi = 0{mb) sums the large 
logarithm a^ \n{Mw/mh) to all orders in perturbation theory. The operator product expansion lead- 
ing to (|]) and the RG improvement amount to a simultaneous expansion in ml/M^, as{Mw) and 
as{mi,) of the b — > ccs amplitude. Tie// in (0) reproduces the leading term in the power expansion 
in ml /M^. 

The second step to predict ATsm involves an operator product expansion at the scale m(,. The 
corresponding formalism has been formulated long ago by the hosts of this conference [|10|]. The 
starting point for the calculation of the widths Th of some 6-flavored hadron H is the optical theorem, 
which relates rj:^ to the absorptive part of the forward scattering amplitude of H. Neglecting CP 
violation in the decay amplitude the optical theorem implies 

r^ oc Jm{H\i fd'^xTneff{x)neff{0)\H). (10) 

Now the Heavy Quark Expansion (HQE) [ [T0| ] is an operator product expansion of the forward scat- 
tering amplitude in ([T^). Schematically 



r^ oc GlY.^f'Pjif'/^b) iH\0,{fi)\H) (11) 



OiA'^co 



Here new Wilson coefficients Fj have appeared. They contain the physics associated with scales 
above the matching scale /i = 0{mi,), at which the HQE is performed. The C/s are local operators 
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Figure 2: Leading-order diagrams for T 
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with dimension dj > 3. The HQE is a simultaneous expansion of Th in Agco/'mb and ^^(mfe). 
Increasing powers of Agc-^/mb correspond to increasing dimensions dj of the local operators Oj. 
To calculate AF from the HQE one must extend the above formalism to the two state system 

{Bs,Bs): 
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The corresponding leading-order diagrams are shown in Fig. 0. ( [T2| ) is matched to local operators 
in analogy to ([TT]): 
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The HQE for the AB = 2 transition in ( [T^ requires four-quark operators involving both the 6-quark 
and the s-quark field, i.e. operators with dimension six or higher. The two dimension-6 operators 
appearing in (|T3]) are 



Q = {biSi)v-A{bjSj)v~ 



Qs = {hiSi)s~p{hySj) 



j)S-P- 



(14) 



In the leading order of QCD the RHS of ( [T3t ) is pictorially obtained by simply shrinking the (c, c) 
loop in Fig. |^ to a point. The Wilson coefficients F and Fs also depend on the charm quark mass rric, 
which is formally treated as a hard scale of order m^, since rric ^ Aqcd- Strictly speaking, the HQE 



in ( [131 ) is an expansion in Aqcd/ yf^l ~ 4m^. For the calculation of F and Fs it is crucial that these 
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coefficients do not depend on the infrared structure of the process. In particular they are independent 



of the QCD binding forces in the external Bg and Bg states in (|T3|), so that they can be calculated in 
perturbation theory at the parton level. The non-perturbative long-distance QCD effects completely 
reside in the hadronic matrix elements of Q and Qs- 

The third and final step in the prediction of ATsm is the calculation of the hadronic matrix 
elements with non-perturbative methods such as lattice gauge theory. It is customary to parametrize 
these matrix elements as 

(B,\Q\B,) = IflXiB 



(B.mB.) = --/|M|.j-^B„ (15) 

In the so called vacuum insertion approximation B and Bs are equal to 1. 



3 Next-to-leading order QCD corrections to AF 



The discussion of AT in sect. ^ has been restricted to the leading order (LO) of QCD [[TTl]. The 
only QCD effects included in this order are the leading logarithms a" \n"'{Mw/mb), n = 0,1,2, .. ., 
contained in the C/s of the effective AB = 1 hamiltonian in (Q). To predict AT with next-to- 
leading order (NLO) accuracy one must first include the corrections of order a""*"^ ln"(M^y/mb), 
n = 0, 1, 2, . . ., to these coefficients [jT^]. Second corrections of order as{mh) must be included in 
F and Fs [|13]]. This step requires the inclusion of hard gluon exchange on both sides of ([131). The 
corresponding diagrams are depicted in Fig. |^. The motivations for this cumbersome calculation 
are 

1 . to verify the infrared safety of F and Fs, 

2. to allow for an experimental test of the HQE, 

3. a meaningful use of lattice results for hadronic matrix elements, 

4. to reduce the sizable yU-dependence of the LO, 

5. a consistent use of A^, 

6. the large size of QCD corrections, typically of order 30%. 



3 Next-to-leading order QCD corrections to AF 
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Figure 3: QCD corrections to AF. NLO corrections to diagrams with the penguin operators Q3 
are negligible. 



We will now comment on these points: A necessary condition for the validity of any operator product 
expansion is the disappearance of all infrared regulators from the Wilson coefficients. From our 
explicit calculation we have verified that this is indeed the case at order a^. We found IR- singularities 
to cancel via two mechanisms: 

• Bloch-Nordsiek cancellations among different cuts of the same diagram, 

• factorization of IR-singularities, which end up in {Bs\0\Bs), {Bs\Os\Bs). 

Early critics of the HQE had found power- like infrared divergences in individual cuts of diagrams. 
In response the cancellation of these divergences has been shown [[T4|], long ago before we have 
performed the full NLO calculation. The second type of IR-cancellations occurs between the dia- 



grams in the first and second row of Fig. |3|. Thus when the external meson states in ( [T3| ) are replaced 
by quark states, both sides of the equation are infrared divergent. Yet the IR-divergences factorize 
rendering F and Fs infrared safe. Point 2 above addresses the conceptual basis of the HQE, which 
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is sometimes termed quark-hadron duality. It is not clear, whether the HQE reproduces all QCD ef- 
fects completely. Exponential terms like exp(— Kmb/AQCD)^ for example, cannot be reproduced by 
a power series. [ |T3| ] The relevance of such corrections to the HQE can at present only be addresses 
experimentally, by confronting HQE-based predictions with data. The only QCD informations con- 
tained in the LO prediction for AF are the coefficients of a" In" My/, associated with hard gluon 
exchange along the Pl^-mediated h — > ccs amplitude. The question of quark-hadron duality, however, 
has nothing to do with these logarithmic terms. A meaningful test of this aspect of the HQE therefore 
requires a NLO calculation, which includes non-logarithmic terms of order a^. In view of the suc- 
cess of the HQE in accurately measured B physics observables it is conceivable that the uncertainty 
due to violations of quark-hadron duality is well below the uncertainty from the non-perturbative 
calculation of the hadronic i? -parameters. At present lattice calculation of B and Bg [|l6l] are only 
possible in the quenched approximation, neglecting the effect of dynamical fermions. Unquenched 
calculations of Jb^ are now available, but still a new subject in the field [[TTt]. The third point in 
our list above refers to the fact that QCD predictions obtained on the lattice must be matched to the 
continuum. This involves the calculation of the diagrams in the first row on Fig. |3| in lattice perturba- 
tion theory. A meaningful prediction for AF with a proper cancellation of the renormalization scale 
and scheme dependences between F,Fs and B,Bs then requires a full NLO calculation. The renor- 
malization scale fi is an unphysical parameter and observables do not depend on /i. The truncation 
of the perturbation series, however, introduces a ^u-dependence, which diminishes order-by-order in 
Ug- As mentioned in point 4, the LO result for AF suffers from a huge scale dependence, which is 
substantially reduced in the NLO prediction. Further a LO calculation cannot use the fundamental 



QCD scale parameter A^ QTSO, which is an intrinsic NLO quantity. Finally, as mentioned in point 
6, the calculated QCD corrections are sizable, of the order of 30%, and therefore necessary to keep 
up with the precision of the forthcoming experiments. 

Including corrections of order Aqcz)/"^6 [ |T9| ] to ( 1131 ) we predict [ [l3l ] 

with B and Bs defined in the MS -scheme at ft = rrih. With [jT^] 

B{ii = mb) = 0.80 ±0.15, Bs{fi = mb) = L19±0.20 



one finds 



V F Jbs V245MeV 



2 



(0.162 ± 0.041 ±???). (16) 



The questions marks address the unknown error from the quenching approximation. 
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If Ar is found below the Standard Model prediction, it will be interesting to find out, whether 
this is due to a breakdown of the HQE prediction in ( [T^ ) or a new CP- violating phase in (Q). To 
this end we note that one can determine (p without using the theory prediction for AF^a/, even from 
untagged data alone [H,|2^]. Further the HQE prediction for other width differences, e.g. between 
the B^ and B^ or between the Bg and B^ mesons, involve a similar structure than the prediction 
for AF. The corresponding diagrams are similar to those in Fig. || and Fig. |], but involve Ai? = 
transitions [[l0|, [191]. The width difference between B^ and Bd is insensitive to new physics and 
therefore directly tests the HQE and the lattice calculations of hadronic matrix elements. The small 
width difference between Bs and B^ is mildly sensitive to new physics from penguin contributions 
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